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ABSTRACT 


The  theory  for  the  flow  of  a  weakly  ionized  gas  through  a  parallel  - 
plate,  continuum,  electrostatic  probe  is  developed.  The  flow  is 
separated  into  three  distinct  regions:  (a)  the  inviscid,  neutral  core 
where  electron  conduction  maintains  the  continuity  of  current  between 
the  two  plates;  (b)  the  viscous,  quasi-neutral  boundary  layer  in  which 
the  charged  particle  flow  is  similar  to  ambipolar  diffusion;  and  (c)  the 
one -dimensional,  collision  dominated,  space-charge  sheath.  Analyti¬ 
cal  solutions,  matched  at  the  boundary  of  each  region,  are  presented 
for  the  electron  temperature  in  equilibrium  with  the  gas  temperature 
and  for  the  electron  temperature  constant  at  its  free-stream  value. 

A  criterion  is  given  which  may  be  used  to  determine  whether  electron 
thermal  equilibrium  exists  through  the  boundary  layer.  It  is  shown 
that  the  sheath  voltage  drop  comprises  approximately  60  percent  of  the 
total  plate  voltage  drop.  The  results  also  show  a  very  well  defined 
saturation  current  for  the  double  probe  and  that  this  current  is  con¬ 
trolled  by  ion  diffusion  through  the  boundary  layer.  Expressions  are 
developed  from  the  solutions  which  allow  the  use  of  experimental  data 
to  determine  the  free-stream  electron  density  and  temperature. 
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NOMENCLATURE 

a  Defined  by  Eq.  (67) 

A  Defined  in  Eq.  (13) 

C  Mass  fraction 

Cp  Specific  heat  at  constant  pressure 

D  Diffusion  coefficient 

E  Electric  field 

e  Electronic  charge 
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f'  Blasius  function,  u/u^ 

J  Defined  by  Eq.  (54) 

Jj  Defined  in  Eq.  (13) 

j  Current  density 

K  Mobility  coefficient 

k  Defined  in  Eq.  (13) 

k  Boltzmann's  constant 

kn  e  Thermal  conductivity 

L  Plate  length 

SL  Width  of  plates 

M  Particle  mass 

m  Ratio  of  mass  fraction,  C/Cro 

M0  Free- stream  Mach  number 

n  Number  density 

Pr  Prandtl  number,  ju  Cp/kn 

p  Pressure 

Q  Effective  collision  cross  section  for  momentum  transfer 

R  Defined  in  Eq.  (13) 

5  Schmidt  number,  n/pD in 

T  Temperature 

u,  v  Velocity  components 

V  Electric  potential 

x,  y  Cartesian  coordinates 

a  Defined  in  Eq.  (13) 

■y  Ratio  of  specific  heats,  coefficients  in  Eq.  (24) 

6  Momentum -energy  loss  parameter 

Boundary-layer  thickness 

A  Defined  in  Eq.  (13) 

Permittivity  of  vacuum 
£  Defined  in  Eq.  (13) 

r I  Transformed  boundary-layer  coordinate  (defined  on  page  20) 
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6  T/Tm 

X  Electron-neutral  mean  free  path 

p  Dynamic  viscosity 

v  Electron-neutral  collision  frequency,  kinematic  viscosity 

P  Mass  density 

<P  Intersection  of  the  slope  of  the  probe  voltage-current 

curve  with  Js  (see  Fig.  11) 

u  Te/T 

£2  Defined  in  Eq.  (47) 

SUBSCRIPTS 

a  Standard  atmospheric  conditions 

e  Electrons 

i  Ions 

n  Neutrals 

o  Edge  of  sheath 

s  Saturation 

t  Total 

w  Wall 

o°  Free -stream  conditions 

Note:  MKS  rationalized  units  are  used  in  this  report. 
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SECTION  I 
INTRODUCTION 


In  recent  years  there  has  been  a  great  deal  of  interest  concerning 
the  use  of  electrostatic  probes  in  plasma  diagnostics.  Much  of  this 
interest  has  been  concentrated  on  extending  the  theory  and  use  of  con¬ 
ventional  Langmuir  probes  (free-molecular  probes)  to  flowing  plasmas. 
Some  of  the  techniques  and  theories  for  these  probes  may  be  found  in 
Refs.  1  and  2.  These  probes  are  practical  and  yield  fair  results  when 
used  in  a  low  density  flow.  As  the  density  increases  the  space- charge 
sheath  becomes  collision  dominated  and  a  continuum  theory  is  required. 

The  term  "continuum  electrostatic  probe"  will  be  used  for  a  probe 
operating  in  a  flow  regime  where  the  sheath  thickness  is  much  larger 
than  an  electron-neutral  mean  free  path.  The  analysis  of  such  a  probe 
must  then  include  the  collision  dominated  sheath  equations. 

Numerical  solutions  have  been  obtained  for  the  nonflowing,  weakly 
ionized,  collision  dominated  space-charge  sheath  for  a  spherical 
geometry  and  for  a  Couette  flow.  Notably,  the  numerical  solution  for 
the  spherical  probe  was  obtained  by  Cohen  (Ref.  3)  and  Radbill  (Ref.  4). 
The  plane  sheath  for  a  collision  dominated,  weakly  ionized  gas  has  been 
solved  by  Chung  (Ref.  5)  for  a  Couette  flow.  The  solution  by  Cohen  is 
presented  as  an  asymptotic  theory  for  a  spherical  probe  in  the  limits  of 
(1)  probe  radius  to  Debye  length  ratio  large,  and  (2)  T^/Tg— -  0  for 
arbitrary  probe  voltage.  After  reducing  the  pertinent  equations  for  the 
above  limits,  Cohen  numerically  integrated  the  resulting  differential 
equation.  It  should  be  pointed  out  that  in  Cohen's  work  the  sheath  thick¬ 
ness  is  much  larger  than  the  electron-neutral  mean  free  path.  Radbill, 
through  a  different  numerical  technique,  extended  the  solution  for  the 
spherical  probe  to  include  arbitrary  probe  radius  to  Debye  length  ratio 
and  arbitrary  potentials.  The  solutions  of  Cohen,  Radbill,  and  Su  and 
Lam  (Refs.  3,  4,  and  6)  all  agree  in  the  regions  where  the  parameters 
of  interest  are  equal.  One  feature  that  is  worthy  of  note  in  these  solu¬ 
tions  is  the  failure  to  obtain  a  saturation  current.  A  possible  explanation 
of  this  behavior  is  the  penetration  of  the  electric  field  into  regions  far 
away  from  the  probe.  Since  the  analyses  have  been  for  a  nonflowing 
plasma  the  field  penetrates  farther  with  increasing  probe  potential  in 
order  to  maintain  continuity  and  there  is  no  sharply  defined  space -charge 
sheath  edge.  The  solution  of  Cohen  served  as  a  basis  for  an  analysis  of 
a  flowing  plasma  over  an  arbitrary  body  by  Lam  (Ref.  7).  In  his  analysis 
Lam  investigated  the  probe  characteristics  for  a  model  very  similar  to 
that  Chung  used  in  his  analysis  of  the  Couette  and  stagnation  point  flow. 
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Although  Lam's  solution  is  for  an  arbitrary  body  it  has  to  be  restricted 
to  a  three-dimensional  geometry  for  the  limiting  case  of  zero  velocity 
since  steady-state  two-dimensional  solutions  cannot  be  obtained  unless 
charge  is  supplied  from  some  source.  Both  Lam's  and  Chung’s  solutions 
show  that  the  electric  field  penetrates  far  into  the  flow,  and  it  is  neces¬ 
sary  to  include  a  region  of  near  ambipolar  diffusion  which  is  matched  to 
the  nonconvective  sheath.  Lam’s  analysis  was  for  an  incompressible, 
isothermal  plasma  with  constant  properties  while  Chung  and  his  co¬ 
workers  have  included  compressibility  (Refs.  8  and  5)  and  possible  elec¬ 
tron  thermal  nonequilibrium  (Refs.  9  and  10).  Other  pertinent  analyses 
for  the  stagnation  point  probe  may  be  found  in  Refs.  1,  11,  and  12  which 
are  all  similar  to  Lam's  and  Chung’s  work  in  their  method  of  analysis. 

In  this  paper  the  analysis  is  of  a  double,  parallel-plate  probe  with 
the  method  of  solution  following  that  of  Chung  and  Blankenship  (Ref.  8). 
Analytical  solutions  have  been  obtained  for  the  complete  problem  and  the 
results  verify  the  numerical  work  of  Ref.  8.  Justification  is  given  for 
obtaining  both  an  electron-neutral  thermal  equilibrium  and  frozen  electron 
temperature  solution.  In  Ref.  8  these  are  obtained;  however,  an  incom¬ 
plete  form  of  the  electron  energy  equation  was  used,  and  justification  for 
the  frozen  electron  temperature  solution  is  not  correct,  An  approximate 
relation  is  given  here  for  determining  electron  thermal  equilibrium  by 
considering  a  balance  of  collision  losses  and  the  thermal  conduction  of  the 
electrons.  Relations  are  given  which  enable  the  results  to  be  used  with 
experimental  data  to  yield  information  concerning  the  free-stream  elec¬ 
tron  number  density  and  temperature.  Of  particular  significance  is  the 
approximate  analytical  solution  of  the  collision  dominated,  plane  sheath 
equations  for  arbitrary  values  of  the  parameters  of  interest. 


SECTION  II 

FORMULATION  AND  BASIC  ASSUMPTIONS 


Many  of  the  electrostatic  probe  theories  developed  have  been  for 
spherical  geometries.  Since  it  has  been  shown  by  Lam  (Ref.  7)  that  the 
electric  field  penetrates  far  from  the  probe,  spherical  symmetry  was 
necessary  to  attain  undisturbed  conditions  far  from  the  probe.  Similar 
to  the  spherical  probes  is  the  stagnation  point  probe  proposed  by  Tabbot 
(Ref.  11).  The  inherent  difficulty  in  all  of  these  probes  is  obtaining  free- 
stream  properties  ahead  of  the  bow  shock  from  the  probe  data. 

Because  of  the  problems  discussed  above,  the  author  initiated  a 
study  of  a  double  probe  consisting  of  two  parallel  plates.  Independent  of 
this  Chung  and  Blankenship  (Ref.  8)  published  a  numerical  solution  of  the 
same  geometry. 
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This  type  of  probe  offers  several  advantages.  The  parallel- plate 
geometry  provides  two  very  distinct  electrodes  for  a  very  definite  de¬ 
scription  of  the  voltage  drops.  With  the  plates  at  close  proximity  to 
each  other,  the  electric  field  is  nearly  normal  to  the  surface  every¬ 
where  and  Poisson's  equation  is  essentially  one -dimensional  in  nature. 
The  plates  may  be  diverged  slightly  to  account  for  the  boundary  layer 
buildup  so  that  the  strong  shock  problem  may  be  either  eliminated  or 
reduced  greatly. 


The  formulation  will  be  based  on  the  following  assumptions: 

1.  The  ordinary  viscous  boundary  layer  will  be  analyzed  with  a 
zero  pressure  gradient. 

2.  The  charged  particle  number  densities  are  nearly  equal  in  the 
inviscid  core  and  outer  portions  of  boundary  layer. 

3.  The  sheath  (where  charge  separation  occurs)  thickness  is  much 
smaller  than  the  boundary-layer  thickness  and  much  larger 
than  the  electron-neutral  mean  free  path. 

4.  The  flow  is  assumed  to  be  frozen  in  ionization  and  recombina¬ 
tion.  (Calculations  are  given  which  show  the  approximate 
range  of  validity  for  this  assumption.  ) 

5.  The  ions  are  in  thermal  equilibrium  with  the  neutrals  in  all 
cases  <Tj  =  T)  and  K{/Ke  is  assumed  constant. 

Within  the  framework  of  the  above  assumptions  the  neutral  gas  equations 
become: 

Overall  continuity: 


5pu  5p v 
bx  +  by 


(1) 


Overall  momentum: 


pu 


du 

dx 


+  p  V 


du 

by 


dy(  M 


du  » 
dy  ) 


(2) 


Neutral  gas  and  ion  energy  equation: 


r>  dT  ,  „  _  o  r  o  /  ,  or v  ,  /  ou  \ 

pu^p  dx  P  p  iy  3y  (  dy)+  ^  (  By  ' 

The  conservation  equations  for  the  charge  particles  become: 
Conservation  of  ions: 

bC. 


b  T 


bT^ 


du 


pu 


bC. 

i 

dx 


+  p  v 


dy  Ty 


[pDl  -  PK.C.E  j 


(3) 


(4) 
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Conservation  of  electrons: 


3C 


ac 


pu  ^7^  +  pv 


3x 


ay 


a 

Sy 


j^pD 


T_  3_  /  CeTe 
e  Te  9y  '  T 


)+  PKeCeE] 


Electron  energy: 

aTe 

puCp  Ce7^  + 
e  3x 


ovC  -  pD  -L-  ( 

■  e  M  e  Sy  V 


C  T 
e  e 


<5) 


)  "  PK, 


C  E 
e  e 


!!e 

*y 


+  f  5  velT  nek  (Te  “  T) 


Poisson's  equation: 


3E  e  ,  v 

—  =  —  (n.  -  n  ) 

3y  e  i  e 

o 


ep 

e 


C. 

l 

M7 


(6) 


(7) 


The  analysis  is  broken  down  into  three  regions  (Fig.  1).  The  first 
is  the  inviscid  neutral  core  in  which  the  electric  field  is  constant  and 
the  charged  particles  are  delivered  to  the  boundary  layer  because  of 
their  mobility  and  convection.  The  second  region  is  the  outer  portion 
of  the  viscous  boundary  layer  where  n^  -  ne  and  the  flow  is  characterized 
by  diffusion  similar  to  ambipolar  diffusion.  The  third  region  is  the 
space -charge  sheath  in  which  nj  i  ne  and  the  convection  can  be  neglected 
in  the  equations  of  motion. 


SECTION  III 

SOLUTION  TO  THE  PLANE  SHEATH  EQUATIONS 

Analysis  of  Eqs.  (4)  through  (7)  is  very  difficult  when  ni  i  ne,  but 
the  complexity  of  the  equations  is  greatly  reduced  if  the  convection  terms 
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are  neglected.  In  the  outer  portion  of  the  sheath  where  quasi-neutrality 
is  attained,  the  convection  does  not  contribute  to  the  net  current  to  the 
wall.  The  region  where  charge  separation  occurs  is  very  much  thinner 
than  the  viscous  boundary  layer,  and  convection  becomes  negligible. 
Therefore  the  convection  terms  will  be  neglected  in  the  sheath,  and  the 
equations  reduce  to  those  of  a  plane,  collision  dominated  space-charge 
sheath. 


As  pointed  out  in  Ref.  9,  thermal  nonequilibrium  may  exist  between 
the  electrons  and  neutrals  in  the  boundary  layer  and  sheath.  It  is  shown 
in  Section  IV  under  what  conditions  this  might  occur  in  argon  and  air. 
There  are  two  limiting  conditions  which  cover  a  wide  range  of  actual 
flows:  equilibrium  electron  temperature  (Te  =  T)  and  frozen  electron 
temperature  (Te  -  T6tj;).  Hence,  the  sheath  equations  will  be  solved  for 
these  two  extremes.  Under  these  conditions  the  electron  energy  equa¬ 
tion  becomes  extraneous,  and  the  plane  sheath  equations  are: 


3_ 

9y 


ac. 

_ i 

3y 


o  K.C.E 

l  l 


dy 


dE 

Sy 


n  T  5  /  e  e\ 

I  p  De  T0  9y  V  T  ) 

[C.  Cl 
l  e 

m  .  r 

1  e  _ 


+  p  KCE 
e  e 


0 


The  first  two  of  these  equations  may  be  integrated  to  give: 


e 


Lp 


p  K.C.E 

li 


(8) 

(9) 


(10) 


(11) 


e 

p  D  X-  1 

C  T 

f  %  e  ' 

)  +  p  K  C  E 

M 

e 

e  Te  8y  1 

L  T  , 

1  e  e 

=  J, 


(12) 


where  both  currents  are  taken  as  positive  for  particle  drift  toward  the 
wall.  Since  the  equations  are  identical  to  those  analyzed  numerically 
in  Ref.  8,  the  same  notation  will  be  used  and  the  following  parameters 
defined: 
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c. 


a . 


i  ,e 


C  . 


i  >e. 


«  -  s 


y 

1  f  EL. 

a  J  po 


dy 


j  .Ap 
l  Ko 


J,  ~ 

1  e  n 


pu- 


-/ 


dy 


„  e  Ki 


A  = 


2  n  2 
e  e  A 
o 


e  o  ^  To 


R  = 


eo  E 


e  ru  A 
eo 


(13) 


Substitution  of  Eq,  (13)  into  Eqs.  (10)  through  (12)  gives 


da. 


dC 


-  ARai  = 


dC 


(a  to)  +  ARa  =  J,  k 
e  el 


(14) 

(15) 


dR 

dC 


=  a .  -  a 
l  e 


(16) 


In  normalizing  the  above  equations,  it  has  been  assumed  that 
ni^  =  neQ.  Equation  (16)  shows  that  this  is  not  exactly  true  (oq^  /  ) 

since  the  electric  field  is  decaying  through  the  sheath.  However,  to  the 
first  approximation  (quasi-neutral  approximation)  this  is  correct  so  that 
the  boundary  conditions  at  ?  =  1.  0  are  taken  as 

a  =  a.  =1.0 

e  i 

o  o 

The  boundary  conditions  at  the  wall  (£  =0),  although  not  exact,  are 
taken  as 

a  =  a .  =0 


A  discussion  of  these  boundary  conditions  is  given  in  Section  3.1. 
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Equations  (14)  and  (15)  can  be  integrated  formally  to  give  oq  and  ae 
as  functions  of  the  normalized  electric  field  R  and  the  independent 
variable  £  . 


(17) 


The  quantity  contains  the  factor  S/pq  and  is  weakly  dependent  upon 
the  temperature.  Since  the  main  contribution  to  the  integrals  is  near 
the  wall  and  the  temperature  does  not  vary  greatly  across  the  thin  sheath, 
the  ion  Schmidt  number  is  assumed  constant  through  the  sheath  and  equal 
to  its  value  at  the  wall.  This  means  that  Ji  will  be  taken  as  constant 
through  the  sheath  and  equal  to  its  wall  value. 


If  the  correct  relationship  between  R  and  ?  were  known,  the  inte¬ 
gration  of  Eqs.  (17)  and  (18)  could  be  carried  out  to  yield  the  solution 
since  u  is  a  known  function  of  i  .  The  exact  form  of  R  requires  the 
solution  of  the  equations  to  be  known;  however,  a  reasonable  approxima¬ 
tion  to  R  will  give  an  approximation  to  oq  and  ae.  These  expressions  for 
cq  and  ae  in  turn  could  be  substituted  into  Eq.  (16)  and  integration  would 
provide  an  improved  approximation  to  R.  This  procedure  is  simply  the 
method  of  successive  approximations  and  is  commonly  used  in  solving 
nonlinear  differential  equations.  The  practical  success  of  this  method 
depends  almost  entirely  upon  being  able  to  obtain  a  very  good  first 
approximation  to  the  function  in  question.  The  choice  of  the  first  approxi¬ 
mation  will  be  considered  when  treating  specific  cases. 

Two  cases  will  be  solved  which  are  of  interest:  (1)  equilibrium 
electron  temperature,  and  (2)  constant  electron  temperature. 
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3.1  EQUILIBRIUM  ELECTRON  TEMPERATURE 


The  condition  of  the  electron  temperature  equal  to  the  ion  tempera¬ 
ture  gives  (j  =  1.0.  Equations  (14)  and  (15)  may  then  be  added  and 
combined  with  Eq.  (16)  to  give 

3C  (ai  +  “e5  "  AR  HT  =  J1  (1  +  (19) 

Integration  of  Eq.  (19)  yields: 

(a-  +_ae)  -  2  -  ~  (R2  -  Rq2)  =  (1  +  k) <£  -  1)  (2o) 

where  the  subscript  o  denotes  conditions  at  ?  =  1.  0.  Equation  (20)  is 
an  exact  algebraic  equation  between  all  the  variables  of  the  problem. 


Subtracting  Eq.  (14)  and  (15)  and  using  Eq.  (16)  gives: 


^  -  AR  (a.  +  a  )  =  J  (1  -  k)  (21) 

dC2  1 

This  equation  along  with  Eq.  (20)  represents  a  second-order,  nonlinear 
differential  equation  for  the  electric  field  R.  An  approximation  to  R  is 
obtained  by  neglecting  d2R/d?2,  which  from  Poisson's  equation  implies 


da 
_ e 

dC 


da. 

_ l 

dC 


(22) 


In  the  outer  portion  of  the  sheath  where  the  flow  is  quasi-neutral, 
oq  =  ue  and  Eq.  (22)  is  a  good  approximation.  As  the  wall  is  approached 
the  electric  field  becomes  stronger  and  the  ions  become  mobility  limited 
such  that  doq/d£  becomes  very  small.  The  electrons  are  being  repelled 
by  the  electric  field  and  their  density  becomes  low.  From  Eq.  (15)  we 
have 

da 

— -  =»  -  AR  a  “0 
e 


for  a  =“  0. 
e 

Therefore,  the  second  derivative  of  R  is  small  over  a  large  portion  of 
the  sheath. 

Making  the  approximation  of  Eq.  (22)  in  Eq.  (21)  gives: 

J  i  - 

(a  +  a  )  - - ±  (1  -  k)  (23) 

S  AR 
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It  should  be  noted  that  Eq.  (2  3)  at  £  =  1.0  gives  the  same  quasi-neutral 
field  as  that  obtained  by  letting  oq  =  in  the  outer  region  of  the  sheath. 
This  field  is  discussed  by  Chung  (Ref.  5)  as  becoming  asymptotically 
correct  for  large  A.  Substitution  of  Eq.  (23)  into  Eq.  (20)  gives: 

C  “  r-L  -  r2  s'1  -  r 3  R2  (24) 

where 

2  +  ARq2 

(1  +  k)  2 J x  (l  +  k) 


^2 


1  -  E 

1  4  k 


^3 


A 

2J1  (1  4  k) 


Equation  (24)  is  an  approximate  expression  for  ?  in  terms  of  R.  A 
comparison  of  Eq.  (24)  with  a  numerical  solution  from  Ref.  5  is  given 
in  Fig.  2,  where  Jj  and  k  were  taken  from  Ref.  5,  and  shows  that  the 
approximation  to  R  is  in  fact  quite  good.  The  highest  order  derivative 
of  R  has  been  neglected  and  two  constants  of  integration  have  been  lost 
The  result  is  that  Eq.  (24)  does  not  satisfy  all  the  necessary  boundary 
conditions. 


There  is  an  important  feature  of  this  solution  which  will  be  utilized 
in  the  constant  electron  temperature  solution  and  can  be  discussed  at 
this  point.  Since  the  electric  field  at  the  probe  surface  is  large  com¬ 
pared  to  its  value  at  the  sheath  edge,  the  ions  tend  to  become  mobility 
limited  near  the  surface.  This  means  that  physically  it  is  the  drag  force 
between  the  ions  and  neutrals  that  retards  the  ion  motion  and  that  the 
concentration  gradients  are  no  longer  important  in  determining  the  par¬ 
ticle  flux.  Equation  (14)  shows  that  for  this  condition  we  have 


a. 


w 


AR 


(25) 


w 


for  the  surface  boundary  condition  {£  =  0).  The  approximation  in  Eq.  (23) 
gives  at  £  =  0 


a .  4  a 

ew 


AR 

w 


(1  -  k) 


(26) 
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which,  for  the  range  of  k  of  interest,  is  very  nearly  equal  to  Eq.  (25). 

It  was  pointed  out  previously  that  Eq.  (23)  becomes  asymptotically 

correct  at  £  -1.0  for  large  A;  then  in  view  of  the  above  discussion, 

the  approximation  in  Eq.  (23)  can  be  expected  to  be  quite  good  through 

the  entire  sheath,  and  making  cq  =0  causes  only  a  small  perturbation 

w 

in  the  electric  field  near  the  probe  surface. 


By  substituting  Eq.  (24}  into  Eqs.  (17)  and  (18),  with  u  =  1.0,  and 
integrating  we  have 


a . 

l 


1/3 


(27) 


a  = 
e 


-  J^k  t 


r2k 


2h 

3A2 


) 


(28) 


where 


and 


t 


2^AR3 


3 
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Equations  (27)  and  (28)  may  now  be  made  to  satisfy  all  the  boundary 
conditions  at  £  =  1.  0  and  £  =  0  by  evaluating  the  constants  (Jj  and  k)  in 
the  equations.  Instead  of  specifying  the  current  or  voltage  across  the 
sheath,  the  ratio  of  currents  (k)  wrill  be  chosen  as  the  parameter;  then 
Jl>  Rw,  and  RQ  will  be  found  by  solving  the  system  of  equations  corre¬ 
sponding  to  Eqs.  (27)  and  (28)  evaluated  at  ?  =  1.0  and  Eq.  (24)  evalu¬ 
ated  at  £  =0. 

Equations  (27)  and  (28)  may  be  integrated  in  terms  of  the  general 
Kummer  function;  however,  it  is  possible  in  most  cases  to  make 
approximations  consistent  with  the  values  of  the  parameters  of  interest 
which  permit  these  integrals  to  be  expressed  in  simpler  forms.  Such 
approximations  and  the  corresponding  integrals  are  shown  in  Appendix  II 
(Eqs,  (11-10),  (II- 11),  and  (11-12)). 

3.2  CONSTANT  ELECTRON  TEMPERATURE 

A  constant  electron  temperature  means  that  u  is  not  constant  so  it 
is  not  possible  to  derive  an  approximate  expression  such  as  Eq.  (24). 

Examination  of  Eqs.  (17)  and  (18)  shows  that  the  ion  equation  has 
the  same  form  as  in  the  equilibrium  case  and  the  electron  equation  has 
been  changed  by  the  appearance  of  the  gas  temperature.  It  may  be  ex¬ 
pected  from  the  form  of  these  equations  that  the  functional  form  of  R  and 
of  the  product  R/io  cannot  be  very  much  different  from  that  of  the  equi¬ 
librium  case  as  given  in  Eq.  (24).  Therefore,  the  functional  form  of  £  as 
given  in  Eq.  (24)  is  taken  as  the  first  approximation  relating  ?  to  R, 

C  =  r1  -  r2R~1  -  r3R2  (29) 

The  coefficients  (t's)  in  this  equation  are  not  the  same  as  those  in  Eq.  (24) 
for  the  equilibrium  solution.  These  coefficients  are  found  by  requiring 
that  the  approximation  for  R  agrees  with  the  true  R  at  the  end  points  and 
requiring  the  slope  at  the  wall  to  be  that  given  for  the  ion  mobility  limited 
condition.  These  requirements  yield  the  following  system  of  equations: 

0  *  -  ^V1  -  ^w2 


1-0  =  y  -  y  R  1  -  y  R  2 

J-  *  O  a  O 
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Solving  this  system  of  equations  for  y^,  and  73  yields: 
2 


AR 


r,  = 


w 


3  r2 


2J1(1  +  k) 


2  R 


w 


(30a) 


AR 


r  o  ~ 


w 


2JX (1  +  k)  0 


1  - 


2J  (1  +  k)  R 

x  O 


2  „ 


AR. 


w 


R. 


w  J 


1  -  -  ^ 


3  h,  1  V  n 

!,.h,3 

W  J 


(30b) 


r3 


2J, (1  +  k)  2R  ^ 
1  w 


(30c) 


Since  the  sheath  is  thin,  a  linear  variation  in  the  gas  temperature  will 
be  used. 


-  =  —  ■[  0  +  (0  -  0  )  c"l 

I  w  o  w  bf 

03  03  J 


(31) 


Using  Eqs .  (29)  and  (31)  the  resulting  integration  of  Eq.  (18)  cannot 
be  carried  out  in  a  closed  form.  Let  F  be  defined  by 


1  d£  _  F  6o  ~  * 

co  dR 


(r2  -  r2  R"1  -  r3  R2) (r2  R  2  -  2 yz  r) 


®  w  _o 

+  —  (^2  K  -  ^3  R) 

CO 

An  approximate  expression  is  now  chosen  for  F  such  that  Eq.  (18)  may 
be  integrated.  We  will  use  a  curve  fit  to  F  of  the  form 


F  *  G  =  02R~2  -  283R 


(32) 

Since  the  main  contribution  to  the  integral  in  Eq.  (18)  occurs  near 
£  =  0,  G  will  be  made  to  agree  with  Fat  ?  =0.  Also,  since  F  appears 
inside  the  integral  of  the  exponent,  G  will  be  made  to  agree  with  F  on 
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the  average  over  the  interval  of  integration.  These  conditions  corre¬ 
spond  to 


and 


f]^=0 


GdR 


The  resulting  values  of  the  constants  are: 


g2 


9  AR  R 
w  w  c 

2J,  (1  +  k)  a) 

±  £> 


1  (9o  +  9w)2Jl(1  +  £) 


R. 


29 


w 


AR., 


w 


tv 


3  o 


-|  W  v 


1 

+  — 

R  2 

w 


R 


Tl 


R 


w  J 


A  9 


B3  - 


w 


2«J,  (1  +  k)tu 

1  a 


2R 


w 


Substituting  Eq.  (29)  into  Eq.  (17)  and  Eq.  (32)  into  Eq.  (18)  and  inte¬ 
grating  gives: 


Ar, 


a±  =  -  Jl  t 


2  r. 


1/3 


.*( 

'  3A^  ' 


Vo  A 


3t 


A^„  +  1 


+  1 


e  t  dt 


(33) 


w 
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ui  ae  = 


-  J^kt 


e 


^  ( 


^3 

3A2 


1/3 

) 


where 


2^3A 

3 


R 


3 


(34) 


Another  relation  besides  Eqs.  (33)  and  (34)  is  needed  in  order  to 
determine  the  three  unknowns,  Jj,  R0,  and  Rw.  This  relation  may  be 
derived  in  a  manner  similar  to  that  used  in  obtaining  Eq.  (20).  Adding 
Eq.  (15)  to  Eq.  (14)  gives: 

—  (a  +  a  a>)  -  AR(a .  -  a  )  =  J,  (1  +  k)  (35) 

dC  1  x  e  i 

Substituting  in  Poisson's  equation  and  integrating  gives: 

(at  +  oe  aj)  -  (1  +  coo)  -  |  (R2  -  Rq2)  =  JjCl  +  k)  (C  -  1) (36) 


The  complete  sheath  solution  is  now  obtained  by  letting  k  be  a 
parameter  and  solving  the  three  equations  obtained  by  evaluating  Eqs.  (33) 
and  (34)  at  ?  =1.0  and  Eq.  (36)  at  ?  =  0  for  the  three  unknowns,  J j,  R0, 
and  Rw.  Having  obtained  these  unknowns,  Eqs.  (29),  (33),  and  (34)  can  be 
used  to  compute  distribution  of  the  quantities  across  the  sheath.  Approxi¬ 
mations  to  Eqs.  (33)  and  (34)  are  given  in  Appendix  II  [Eqs.  (11-16), 

(11-17),  and  (11-18)]  which  makes  the  sheath  solution  more  amenable  to 
numerical  calculations. 


3.3  SHEATH  VOLTAGE 


The  sheath  voltage  can  be  found  by  integrating  the  electric  field 


dy 


o 
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which  becomes  for  the  sheath  quantities 

1.0 


kT 


(V  -  V  )  = 
o  w 


9  R  dC 


(37) 


Using  the  linear  variation  for  6  from  Eq.  (31),  the  voltage  becomes 

1.0 

J 

e 


kT 


(vo  -  V  -  -A  « 


°  "  V  / 


C  R  dC 


1.0 

* 

-  A  0_.  /  R  d£ 

0 


’■/ 


(38) 


which  may  be  integrated  using  Eq.  (24)  for  R.  Finally,  we  have  the 
following  expression  for  the  sheath  voltage: 


kT 


(v  -  v  )  =  -A  (e  -  e  ) 

o  w'  o  w 


[^2  ,n-£ 


+  ^2(r-r) 


o  w 


5^2  (R  3  -  a  3)  + 

„  o  w 


^  <C  -  O 


(39) 


2^3 

+  -5-  (R 


o5  -  O  ]  -  A  ew  [  ^2  ^ 

W 


*3  «„3  -  R.3)] 


W 
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where  the  y's  are  given  in  Eq.  (24)  for  the  electron  equilibrium  case  and 
in  Eq.  (30)  for  the  constant  electron  temperature  case. 


SECTION  IV 

BOUNDARY-LAYER  SOLUTION 


The  region  in  the  viscous  boundary  layer  but  outside  the  space - 
charge  sheath  is  quasi-neutral  (n^  =  ne).  Observing  this  it  is  possible 
to  eliminate  the  electric  field  from  Eqs.  (4)  and  (5)  to  give. 


pu 


pv 


3m 

dy 


d_ 

By 


pD.  K.  K  , 
_ 1  1  e  oin 

Ik.  K.  +  K  dy 
L  1  1  e  J 


pDe  Ki  Ke  I  T  ^  TeNl 


where  m  =  and  mj  =  mg.  Since 

k.  m  1/2 

IT  '  M  / 


e  1 

then  Kj  «  Ke  and  Eq.  (40)  becomes. 


dm 

pu  -  4- 

dx 


PV  —  =  —  jpD  —  I"  m  (  1  +  —  )  1  | 
By  By  f  1  Sy  L  T  > 


(40) 


(41) 


which  is  the  convection-diffusion  equation  for  the  charged  particles.  It 
can  be  seen  that  Eq.  (41)  is  similar  to  the  neutral  gas  energy  equation, 
and  similar  solutions  may  be  found  under  certain  assumptions.  Before 
looking  for  similar  solutions  it  is  necessary  to  look  at  the  electron 
energy  equation  since  Eq.  (41)  is  coupled  to  it. 
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4,1  DETERMINATION  OF  ELECTRON  TEMPERATURE 


Equation  (6)  can  be  written  as: 

5  3T  r  T  S  C-T 

2  e  ax  [  e  e  Te 


(-V)  -pw] 


[ 


K  i  ST 
5  k  e 


*1  ~  .  3T  » 

—  E  1  -  —  (  ke  —  ) 

M  j  3v  '  av 


M  Sy  M  J  5y  '  ~  Sy 

0  © 


o  M 

+  “  4  ve  —  nek  (Te  “  T)  <42> 

n  c  C  C 


where  Ref.  13  gives, 


k  = 
e 


15k#)  (r)V2 


4  M  v  T 
e 


Since  the  diffusion  velocities  are  small,  these  terms  are  neglected  and 
Eq.  (42)  reduces  to. 


_  aT  -  st 

5  ,  e  5  ,  e 

—  n  ku  -  +  —  n_kv 


e  3x  2  e 


Sy 


2>t  q 

(k  — sw  1 
V  e  I  2 


Sy  '  Sy 


Me 

6v  —  n  k  (T  -  T) 
e  M  e 


(43) 


The  last  term  in  this  equation  is  similar  to  the  source  term  in  the  con¬ 
servation  of  species  in  reacting  boundary  layers.  In  order  to  get 
similar  solutions  for  Eq.  (41)  it  is  necessary  that  Te  can  be  expressed 
in  a  similar  solution.  It  is  obvious  that  if  the  collision  term  dominates, 
then  Te  =  T  and  Eq.  (41)  has  similar  solutions.  Also,  if  the  collision 
frequency ‘is  small  so  that  the  collision  term  is  negligible,  then  the 
electron  temperature  is  frozen  through  the  boundary  layer.  Similar 
solutions  are  again  possible  for  Eq.  (41)  since  T  can  be  found  from  the 
solution  of  the  neutral  gas  energy  equation.  This  case  is  difficult  to 
solve  in  general  since  the  equation  is  highly  nonlinear.  The  possibility 
of  electron  thermal  nonequilibinum  is  known  and  has  been  investigated 
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for  a  stagnation  point  probe  by  Chung  and  Mullen  (Ref.  10).  Similar 
solutions  might  be  obtained  for  Eq.  (43)  if  some  type  of  variation  in 
the  gas  properties  were  allowed. 

Let  us  assume  a  constant  electron  thermal  conductivity  and  write 
Eq.  (43)  as: 

8T„  0  ..  T_  „  1/2  92Te 

dy2 


-  M 

+  -  —  (T  -  T)  (44) 

5  e  M  e 

In  order  to  determine  the  important  terms  in  Eq.  (44)  we  define. 


DO 


where,  in  the  usual  order  of  magnitude  analysis. 


a_ 

axf 


^  0 


and  6g  is  the  viscous  boundary-layer  thickness.  The  order  of  magnitude 
of  the  terms  in  Eq.  (44)  are  then. 


u' 


ST  ’ 
e 

ax' 


ST 


+  v’ 


ay’ 


3  1  (  M_  ) 

2  R-  '  ' 


1/2  52T  - 
e 

3y’2 


-  T 

+  -  6v  - - —  (  T  ’  -  T  ) 

5  e  u  M  V  e  ' 


(45) 


(1) 


(1)  +  (6  )  -1—  -  )1/2  1 

(6B>  R  ^  '  (6n2) 


e  "e  '~B 

+  collision  term. 
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The  thickness  of  the  viscous  boundary  layer  is 


B 


R 


x 


so  the  electron  thermal  conduction  term  is  of  order 

1/2 


0  (  -) 
V  M  / 


compared  to  unity  for  the  convection  terms.  The  magnitude  of  the  col¬ 
lision  term  depends  on  the  energy  exchange  between  the  electrons  and 
neutrals.  Equation  (45)  is  essentially  a  balance  between  the  thermal 
conduction  and  the  collision  loss  terms  and  can  be  written  as: 


32T 


e  5y2 


e  =  -  -  6v 


M 


M 


n  k  (T  -  T) 
e  e 


(46) 


When  the  right  side  dominates,  Te  =  T;  when  the  right  side  is  negligible, 
Te  -  constant  through  the  boundary  layer.  To  develop  an  approximate 
criterion  for  electron  thermal  equilibrium  we  let 


S2T 


ay' 


-  Te 


w 


(  T  T 
\  e  -  e 

'  DO 


J 


6b2 


25  x2 


R 


x 


Equation  46  then  reduces  to 

— -  -  -  f(i  -  n)  + 

T  2  l 


(i  -  n)  +  4ne 


»]1/2i 


(47) 


where 


0w  " 


w 


SI  = 


10  8v  x6 
_ e  w 

u 


M  3/2 
(  —  )  =  7.15  x  10 

M 


x5Q  P 

21  en  co 


M 


_  0 


W 
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Equation  (47)  is  shown  in  Fig.  3  for  two  representative  wall  tempera¬ 
tures.  The  magnitude  of  ft  depends  on  gas  flow  conditions  and  the  type 
of  gas  being  used.  For  argon  6  =  10  and  Qen  -■  7  x  10"21  ni2j  and  for 
nitrogen  or  air  6  probably  lies  between  10  and  100  and  Qen  s  10"19  m2. 
Using  the  typical  flow  conditions  of  =  4200°K,  =  0.  1  atm;  M„  =1, 

and  0W  =  0.  1  -  0.3  the  corresponding  values  of  ft  for  argon  and  air  are 
shown.  It  is  shown  that  for  argon  the  electron  temperature  is  nearly 
constant  and  for  air  (6  =  100)  the  electron  temperature  is  nearly  in 
thermal  equilibrium  with  the  gas.  Hence,  the  two  extremes  of  equili¬ 
brium  Te  and  constant  Te  are  very  representative  of  many  flows  of 
interest.  It  should  be  pointed  out  again  that  these  expressions  were 
developed  under  the  assumption  that  the  ionization  and  recombination 
are  frozen. 


4.2  SOLUTION  FOR  THE  CHARGED  PARTICLE  CONVECTION-DIFFUSION  EQUATION 


Equation  (41)  can  be  solved  approximately  for  the  two  cases  of  the 
electron  temperature  if  the  quantity  m  is  assumed  to  vary  slowly  along 
the  edge  of  the  sheath.  The  equation  can  be  transformed  by 


T1 


.  U 

-  (  — ) 
'2  v  x  ' 

CO 


1/2  * 
1  / 


—  dy 


t  =  (2u  V  x)1/2  f (n) 


and  letting 


PI* 


P v  _  _  dju  pu  _  df  _  u_ 

p  3x  p  ay  dr|  U 

OO  OO  00 

=  1.0.  Applying  this  transformation  yields 


dn  \S 


1  d  (rng)  ^  t  f  dm  _  Q 


dq  \  S  dil 


dq 


(48) 


where  g  =  (1  +  ub  The  boundary  conditions  on  Eq.  (48)  are 

m  (co)  =1.0 

“  =  mo 

where  r/0  is  the  sheath  thickness. 
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Equation  (48)  can  be  written  in  the  following  form: 

d_  r  1  d  (mg)  1  +  f  d  (mg)  =  mf  din  g  (49) 

dri  L  S  dp  J  g  dp  dp 


For  equilibrium  electron  temperature  the  right  side  of  the  above  equa¬ 
tion  is  identically  zero.  For  constant  electron  temperature  g  becomes 
a  function  of  the  gas  temperature  and  the  right  side  can,  in  general,  be 
neglected.  Neglecting  the  right  side  of  the  above  yields: 

d_fl  dlingn  +  f  dfagl  .  o  (50) 

dp  L  S  dp  J  g  dp 

which  is  to  be  used  as  an  approximation  to  Eq.  (48).  Integration  of 
Eq.  (50)  gives 


mg  -  m0gQ  +  (g^ 


(51) 


The  exponentials  in  the  integrand  of  the  above  equation  have  their  largest 
contribution  near  the  wall  and  decay  as  rj  increases.  The  same  type  of 
behavior  is  seen  in  Appendix  III  in  the  approximate  solution  of  the  Blasius 
equation.  It  is  shown  there  that  reasonable  approximations  can  be  ob¬ 
tained  by  using  the  correct  behavior  of  the  exponentials  near  the  wall  pro¬ 
vided  that  S  and  g  do  not  vary  greatly.  Therefore  the  following  approxi¬ 
mations  are  used  in  Eq.  (51): 

s  -  V  s  * 

and 

f  "  p2 

f  -  - 

2 
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This  gives. 


mg  -  mogo  +  (g_ 


S  S*  /  1 .  f  w  Sw  _  3 

1/ Fs-  i  )-y\s-  ~Bg^-  \  )J 


f  S 


mogo) 


r(k)-r(b%^%3) 


W 


(52) 


The  solution  for  the  equilibrium  electron  temperature  can  be  obtained 
by  letting  g  =  2  everywhere  in  Eq.  (52). 


The  order  of  the  approximation  made  to  Eq.  (48)  by  Eq.  (50)  can  be 
evaluated  by  formally  integrating  Eq.  (49)  and  comparing  the  terms.  It 
can  be  shown  that  if  the  following  integral  is  much  less  than  unity  the 
approximation  is  good: 


I  =  -  — 


mf  dg 
g  dp 


P 

/  M 

*  g 

r\ 


dp 


L  p 


Se 


-/  ft  dll 


dp 


dp 


To  evaluate  I  we  let 


e  -  K  +  (1  -  0W)  f' 


and  from  Eq.  (52) 


m  =3 


2f  ' 

g 


The  integral  then  becomes 


-j-  [ 

J  dp 


dg  1  -  f 


dp 


dp  f 


6g„ 


The  term  neglected  is  on  the  order  of  10  percent  or  less  since  g^  >  2, 
hence  Eq.  (50)  is  a  good  approximation  to  Eq.  (48). 

Equation  (52)  will  be  used  in  the  matching  of  the  boundary  layer  to 
the  sheath  at  the  outer  edge  of  the  sheath. 
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4.3  BOUNDARY-LAYER  VOLTAGE  DROP 


The  electric  field  has  to  be  obtained  from  the  species  equations, 
since  to  the  first  order  -  ne  and  Poisson's  equation  gives  no  informa¬ 
tion  about  the  electric  field.  If  we  observe  from  Eqs.  {4)  and  (5)  that 
the  convective  portions  of  these  are  equal  we  obtain  an  equation  which 
may  be  integrated  from  the  edge  of  the  sheath  outward.  This  can  be  re¬ 
duced  to  the  following  expression  for  the  electric  field: 


E  - 


1 


m 


+ 


p  kT 

00  00 

p  e 


(53) 


where  the  integration  constant  has  been  taken  to  be  J  =  J  (x).  The  J  is 
evaluated  by  considering  conditions  at  the  edge  of  the  sheath  where  the 
effects  of  charge  and  separation  and  convection  are  small;  with  Eqs.  (11) 
and  (12)  J  is  found  to  be 


s  ' 

J  -  (  JS-  )  (j  .  j  )  (54) 

en  V  v  u  /  e  1 

e  co  oo 

OQ 


The  voltage  drop  from  the  edge  of  the  sheath  to  the  edge  of  the 
boundary  layer,  rj  =  p ^ ,  is  given  by  the  integration  of  Eq.  (53), 
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<V5  -  V  * 


-  e) 


e  do 


r  h 
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L  Y| 


6 


9  d^n  m  dp  +  J 
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—  dp 
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'0 

f  e/dfejON  dt) 

J  m  V  dr,  1 


(55) 
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For  the  equilibrium  case  u  =  1  and  Eq.  (55)  becomes 

n. 


kT 


<V6  -  Vo>  - 


/ 


i  —  dri  +  0 
m  dr| 


M 


_  | 

L\m/ 


v  l/2-| 


(56) 


To  evaluate  the  integral  we  shall  use  6  as  given  by  the  compressible 
flow  solution  with  Pr  =  1,  which  is  (Ref.  14), 

e  =  0  +  (i  -  9  )f  +  r  -  1  m  2  f’  (i  -  f') 

w  w  2  “ 


where 


u 


Then  Eq.  (56)  becomes, 
e 


kT 


(V6  -  V  "  -  9w  'tn  % 


+ 


[  -  V  ♦  ^  «.»] 


dm 

dri 


dri 


Y  ~  1 


M 


6  2  /  r  m  1/2  ^ 

!  ^  ^  dr)  +  0  (  T  -i]  ) 

J  m  dq  / 


(57) 


The  largest  contribution  to  the  remaining  integrals  occurs  near  r?  =  rjG 
f  *  dm 

since  — - *- 1  and  decays  exponentially  as  ry—*-v Therefore  any 

approximations  made  in  the  evaluation  of  the  integrals  should  be  very 
good  near  t)0.  Expansion  of  Eq.  (52)  for  q  near  qQ  gives 


m  = 


S  1/3 

=  m  +  (1  -  m  )  f  "  (  —  )  (T|  -  T|  ) 

o  °  w  \  2  '  ° 
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where  we  have  used 


which  is  consistent  with  the  approximate  evaluation  of  Eq.  {52). 
Also, 

f '  =  fj'n  +  o  (ri2) 

If  mQ  «  1  and  t?0  «  1  then 


over  most  of  the  range  of  integration.  The  first  integral  in  Eq.  (57) 
becomes 


•ns 


and  the  second  integral  becomes 


Numerical  integration  of  a  few  cases  has  shown  that  these  approxima¬ 
tions  are  very  good.  Finally,  the  equilibrium  boundary-layer  voltage 
droD  becomes 
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u 

The  nonequilibrium  boundary  layer  with  u  =  in  Eq.  (55)  can  be 
integrated  directly  to  give 
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kT 

00 


CD 

00 


m 


+ 


(59) 


Equations  (58)  and  (59)  will  be  added  to  the  sheath  voltage  drop  to  give 
the  complete  boundary-layer  voltage  drops. 


4.4  INYISCID  CORE  VOLTAGE  DROP 


In  the  inviscid 


region  we  have  from  Eq.  (53) 


kT 

E  =  _ ”  j 


1/2 


(  — ) 
'  2v  x  ' 


(60) 


The  voltage  drop  between  the  two  boundary-layer  edges  is 


where  &  is  the  distance  between  the  plates.  Hence,  across  the  inviscid 
core 


If  the  integration  in  Eq.  (61)  is  from  the  more  negatively  biased  plate 
toward  the  more  positive  plate,  the  voltage  drop  in  Eq.  (61)  will  be 
positive. 


The  voltage  drop  across  the  inviscid  core  may  or  may  not  be  negli¬ 
gible  with  respect  to  the  total  probe  voltage.  The  voltage  drop  between 
the  plates,  neglecting  the  inviscid  core,  (Vwl  -  Vw2),  can  be  written  as 


00 


<P. 
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(62) 
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It  is  shown  in  Section  VI  that 


so  that 
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(63) 


and  the  inviscid  core  voltage  drop  may  be  neglected.  When  this  inequality 
is  not  satisfied,  Eq.  (61)  may  be  used  to  compute  the  inviscid  core  voltage 
drop  and  it  must  be  included  in  the  overall  probe  voltage, 


SECTION  V 

MATCHING  OF  THE  SHEATH  AND  BOUNDARY  LAYER 


Solutions  have  been  given  for  the  sheath  and  boundary  layer;  it  re¬ 
mains  to  match  these  at  the  edge  of  the  sheath  to  complete  the  solution. 

In  matching  the  boundary-layer  solution  to  the  sheath  solution,  three 
conditions  can  be  specified.  These  are  taken  to  be  that  the  electric  field. 
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the  number  density  of  charged  particles,  and  the  electric  current  are 
all  continuous  at  r?  =  n0. 


The  continuity  of  current  has  already  been  assured  by  the  choice  of 
the  integration  constant  in  Eq.  (53).  By  comparing  Eq.  (54)  to  Eq.  (13), 
this  corresponds  to 


J 


(64) 


which  is  the  form  to  be  used  here. 


The  continuity  of  the  number  density  of  charged  particles  has  been 
assured  by  specifying  that  cq  and  are  unity  at  the  edge  of  the  sheath. 

By  comparing  Eq.  (53)  with  the  difference  between  Eqs.  (11)  and  (12), 
it  is  found  that  the  electric  field  is  continuous  provided  the  quantity 
dm  . 

is  continuous  at  17  =  nD.  This  will  now  be  done  for  the  two  cases  con¬ 
sidered. 


For  the  boundary-layer  solution,  we  have  from  Eq.  (52) 


The  sheath  conditions  can  be  expressed  from  Eq.  (35) 

/  \  =  j  (1  +  k)  mo  /  dco  v 

'  dp  •  ^  r|  g  g  \  dp  ' 

o  'o  &o  Bo  1  o 


where  aqQ,  =1.0  and  1  +  wQ  =  gQ  have  been  used.  The  last  terms  of 

the  above  equations  are  identical;  the  matching  of  the  derivatives  at 
p  =  rjQ  then  gives 
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(65) 
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The  nondimensional  current  from  Eq.  (64)  then  becomes 

1/3 

j  -  (1  -  ^ r  r  (!r)  -  li 

1  «.  '  *V  0.  L  \L  J 


(66) 


Introducing  the  parameter 
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into  the  solution  we  have  the  relation 


A  „  2 

“  -  % 
a 


(67) 


Substituting  Eq.  (67)  into  Eq. 

(65)  gives 

r, 3  r  i  -  a  g°  „  l 
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1  s. 

from  which  the  sheath  thickness  can  be  found  for  a  and  k  from  the 
sheath  solutions.  Equation  (68)  contains  a  term  with  gQ  which  is  a 
function  of  r)OJ  however  this  term  is  usually  negligible.  For  equilibrium 
gQ  =  2  and  for  constant  electron  temperature 


g. 


CD 

1  +  — 
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where 
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e  + 
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L  w  Q  00  J  w 


For  a  given  set  of  the  parameters  a.  A,  k,  S,  u*  and  <?w  the  current  and 
voltage  through  the  boundary  layer  are  found. 


SECTION  VI 
RESULTS 


Analytical  solutions  are  given  for  the  flow  of  a  weakly  ionized  gas 
between  two  conducting  plates.  The  solutions  are  given  for  the  particular 
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conditions  of  (1)  equilibrium  electron  temperature,  and  (2)  constant 
electron  temperature  through  the  sheath  and  boundary  layer. 

Typical  sheath  solutions  for  the  charge  particle  density  and  electric 
field  profiles  are  shown  in  Figs.  4  and  5.  Figure  4  also  gives  a  com¬ 
parison  with  an  exact  numerical  solution  for  the  same  set  of  parameters 
and  assumptions  as  taken  from  Ref.  5.  The  agreement  is  very  good. 

This  same  type  of  agreement  was  found  for  the  entire  range  of  param¬ 
eters  given  in  Ref.  5,  indicating  that  the  analytical  solution  is  valid  over 
the  entire  range  of  parameters  of  interest.  Figure  5  also  gives  a  com¬ 
parison  with  an  exact  numerical  solution  of  the  nonequilibrium  sheath  as 
given  in  Ref.  8.  The  analytical  results  are  for  a  constant  electron  tem¬ 
perature  and  a  linear  gas  temperature  profile  while  the  numerical  results 
included  an  approximate  electron  energy  equation  and  solved  for  the 
electron  temperature.  The  agreement  between  the  two  solutions  is  very 
good  despite  the  assumed  constant  electron  temperature  in  the  analytical 
solution. 

The  sheath  voltage  as  calculated  by  Eq.  (39)  for  an  equilibrium  sheath 
is  shown  in  Fig.  6.  An  isothermal  and  a  linear  gas  temperature  profile 
are  shown  in  order  to  give  the  effect  of  compressibility  on  the  sheath 
voltage.  Direct  comparison  with  the  numerical  solution  from  Ref.  5  is 
again  made  for  the  isothermal  sheath,  and  the  agreement  is  very  good. 

The  effect  of  the  density  variation  is  to  decrease  the  voltage  required  to 
produce  a  given  current.  Since  the  electrons  are  in  equilibrium  with  the 
gas,  their  energy  is  decreasing  through  the  sheath.  Therefore,  it  takes 
less  voltage  to  decrease  the  electron  current  and  keep  the  same  net  current. 

Having  demonstrated  the  validity  of  the  sheath  solutions,  they  are 
matched  to  the  boundary -layer  solution  to  give  the  complete  solution.  Typi¬ 
cal  profiles  of  the  complete  solution  are  shown  in  Figs.  Ta  and  b.  It  can 
be  seen  that  the  constant  electron  temperature  sheath  is  thicker  than  the 
equilibrium  sheath.  The  electric  field  decays  very  rapidly  through  the  sheath, 
with  the  value  at  the  edge  of  the  sheath  being  approximately  one-tenth  the 
value  at  the  wall. 

Although  the  electric  field  is  weak  in  the  boundary  layer,  the  voltage 
drop  is  not  negligible  because  of  the  much  larger  distance  involved.  This 
is  illustrated  in  Fig.  8  which' shows  the  sheath  and  boundary-layer  voltage 
drops  for  the  case  of  a  constant  electron  temperature  and  an  equilibrium 
electron  temperature.  The  voltage  drop  across  the  sheath  is  seen  to  be 
approximately  60  percent  of  the  total.  Figure  8  also  shows  a  very  well  de¬ 
fined  saturation  current  for  the  probe.  The  equilibrium  boundary-layer 
voltage  is  a  function  of  the  free-stream  Mach  number  through  the  gas  tem¬ 
perature  dependence.  Equation  (5  8)  shows  that  boundary-layer  voltage 
drop  increases  with  increasing  free-stream  Mach  number. 


30 


AEDC-TR-67-71 


The  sheath  thickness  increases  with  a  decrease  in  the  free-stream 
electron  density  since  the  sheath  thickness  is  characterized  by  (not 
equal  to)  the  Debye  length  which  increases  with  a  decrease  in  electron 
density.  This  variation  of  the  sheath  thickness  with  ne  (n6a}  is  pro¬ 
portional  to  the  parameter  a)  is  shown  in  Fig.  9.  It  is  again  noted  that 
the  constant  electron  temperature  sheath  is  thicker  than  the  equilibrium 
electron  temperature  sheath.  A  discussion  in  Appendix  IV  indicates 
that  the  solution  is  valid  (i.  e.  ,  convection  need  not  be  considered)  for  a 
sheath  thickness  r]0  <  1.4.  Figure  10  shows  the  sheath  thickness  to  be 
directly  proportional  to  the  sheath  voltage. 


The  current -voltage  characteristics  for  the  double  parallel -plate 
probe  may  be  constructed  from  the  total  current -voltage  curves  such 
as  the  ones  shown  in  Fig.  8.  If  the  lower  plate  is  at  a  negative 
potential  with  respect  to  the  plasma,  it  will  draw  an  excess  ion  current. 
The  upper  plate  is  then  less  negative  with  respect  to  the  plasma  and  will 
draw  less  ion  current.  From  the  continuity  of  current  we  have  that 
Jjf  (x)  =  -Ju(x);  hence  the  probe  voltage  is  found  by  subtracting  the  volt¬ 
ages  in  Fig.  8  at  J  and  -J.  Figure  11  shows  these  probe  voltages  for 
several  values  of  the  parameters,  Te  / Tro  and  Sw. 


It  takes  a  much  larger  voltage  to  saturate  the  current  when  Te  /T^ 

is  larger  than  one.  In  most  cases  the  saturation  current  may  be  found 

m0g0  _ 

from  Eq.  (66)  by  neglecting  - .  This  gives,  as  k — >- 0, 
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The  free-stream  electron  density  is  found  from  Eq.  (54) 
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(70) 


Knowledge  of  Te  /Tm  is  required  in  Eq.  (69)  to  determine  Jg  for  use  in 

Eq.  (70).  For  most  cases  this  ratio  will  be  unity;  when  the  ratio  is 
different  than  unity  it  can  be  evaluated  from  the  slope  of  the  current- 
voltage  curve  at  zero  current.  The  saturation  voltage  increases  with 
Te  /Tm,  and  the  slope  of  the  characteristic  curve  is  a  function  of 

<X) 

T6oo /Too  as  was  seen  previously.  Using  a  method  given  by  Chung  and 
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Blankenship  (Ref.  8)  a  correlation  for  Te  /Ta  can  be  found.  A  straight 
line  is  drawn  with  a  slope  equal  to  the  characteristic  curve  slope  at 
J  =  0  from  the  origin  until  it  intersects  Js.  This  voltage  is  denoted  <j)s 
in  Fig.  11.  It  is  found  that  $s  can  be  related  to  Te  /T^  by 


T  2.10 

00 


which  is  very  nearly  equal  to  the  value  given  in  Ref.  8,  even  though  the 
saturation  currents  differ  considerably.  It  should  be  noted  that  the 
voltage  drop  in  the  inviscid  core  has  been  neglected  in  these  figures,  but 
it  may  not  always  be  negligible. 

In  order  to  demonstrate  the  use  of  experimental  data  in  Eqs.  (69) 
and  (71),  by  assuming  a  typical  flow  condition  and  probe  dimensions. 

Fig.  12  shows  the  dimensional  current -voltage  curve  one  could  expect 
to  obtain  experimentally.  The  regions  of  interest  are  marked  on  the 
curve.  Also,  to  aid  in  choosing  instrumentation  for  such  a  probe,  the 
saturation  current  is  shown  as  a  function  of  free-stream  electron  density 
in  Fig.  13.  The  parameter  u  is  a  function  only  of  the  neutral  gas  flow 
variables  and  TQ^/Ta>.  Conservative  estimates  of  the  current  density  may 

be  made  by  letting  Tg^/T^  =  1.  0  in  ct.  The  condition  from  Fig.  12  is 
shown  for  comparison. 

Both  Eqs.  (69)  and  (70)  depend  on  the  ion  Schmidt  number.  This  is 
a  result  of  the  saturation  current  being  dominated  by  ion  diffusion. 
Therefore,  a  knowledge  of  ion-neutral  diffusion  coefficient  is  required. 

The  first  approximation  to  the  binary  diffusion  coefficient  is  given  by 
Demetriades  and  Argyropoulos  (Ref.  15)  as 


where  Qin  is  the  effective  collision  cross  section  for  momentum  transfer. 
This  expression  can  be  written  as 

97  t<3/2  .  2 

D.  =  2  x  10~27  (72) 

P  Qin  sec 
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where  T  is  in  °K;  p'is  in  atm  and  Q^n  is  in  m^.  Not  much  information 
is  available  on  ion-neutral  collision  cross  sections.  It  is  noted  {Ref.  16) 
that  the  collision  cross  sections  are  approximately  three  times  the 
neutral-neutral  collision  cross  sections  for  most  species  of  air.  Also, 
the  magnitude  does  not  vary  too  greatly  from  specie  to  specie,  hence 
we  take 


Q. 

in 


4.36  x  10 


-17  -,-1/2 


(73) 


for  the  effective  collision  cross  section.  The  square  root  dependence 
on  temperature  is  a  result  of  using  inverse  fifth-power  law  interactions 
which  are  characteristic  of  a  simple  polarizable  particle.  Using  Eqs.  (72) 
and  (73)  in  the  definition  of  the  ion  Schmidt  gives. 


s  -  -y- 


PDin 


=  6.44  x  107  ^ 


(74) 


where  p  is  in  Kg/m-sec  and  T  is  in  °K.  Figure  14  shows  the  variation  of 
Eq.  (74)  for  argon  and  nitrogen  as  a  function  of  temperature.  The  vis¬ 
cosity  of  the  neutral  gases  was  calculated  using  viscosity  data  and  a 
Lennard- Jones  potential  (Argon  in  Ref.  17  and  Nitrogen  in  Ref.  18).  If 
a  linear  viscosity  law  is  used,  consistent  with  the  pp  =  constant  approxi¬ 
mation,  the  Schmidt  number  is  a  constant.  Clearly,  the  best  available 
cross  section  or  diffusion  data  should  be  used  in  the  probe  theory  when 
interpreting  experimental  data. 


SECTION  VII 
CONCLUSIONS 


An  electrostatic  probe  has  been  analyzed  for  a  flow  regime  in  which 
the  sheath  is  collision  dominated.  The  probe  consists  of  a  double 
parallel-plate  arrangement  with  the  actual  current-carrying  segments 
being  far  from  the  leading  edge.  The  aerodynamic  boundary  layer  is  in¬ 
cluded  in  the  analysis  along  with  the  continuity  and  energy  equations  of  the 
charged  particles. 

Analytical  solutions  have  been  obtained  for  this  problem.  From 
these  solutions  relations  have  been  obtained  which  allow  the  determina¬ 
tion  of  the  free-stream  electron  density  and  temperature  from  experi¬ 
mental  probe  data  (see  Eqs.  (70)  and  (71)). 

It  is  shown  that  the  flow  consists  of  three  regions:  (1)  the  inviscid 
core  where  the  electric  field  is  very  weak  and  the  current  is  by  electron 
conduction  which  serves  to  maintain  continuity  between  the  plates. 
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(2)  the  viscous  boundary  layer  in  which  the  controlling  mechanism  is 
ion  diffusion  and  is  similar  to  ambipolar  diffusion  except  that  there  is 
a  finite  current  flow,  and  (3)  the  space -charge  sheath  across  which  ion 
diffusion  and  conduction  play  equal  roles.  It  is  found  that  this  latter 
region  contains  very  large  electric  field  gradients  and  hence  a  large 
part  of  the  voltage  drop. 

The  saturation  current  is  dominated  by  ion  diffusion,  and  accurate 
knowledge  of  the  ion-neutral  diffusion  coefficients  is  required  for  the 
particular  gas  being  investigated.  If  measurements  of  n6a)  and  Tg^/T^ 
can  be  obtained  from  other  sources,  the  continum  electrostatic  probe 
could  be  used  in  reverse  and  predict  ion-neutral  diffusion  coefficients. 

The  range  of  free -stream  electron  density  for  which  the  theory  de¬ 
veloped  is  applicable  is  discussed  in  Appendix  IV.  Although  the  validity 
of  several  of  the  assumptions  depends  on  the  particular  flow  (such  as  the 
ratio  of  boundary-layer  thickness  to  sheath  thickness)  it  appears  that  the 
probe  may  be  used  in  the  range  of  lO^/cm^  <  ne  <  10^/cm^.  These 

C<X> 

figures  are  approximate  and  depend  upon  the  particular  flow  conditions 
being  considered. 
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A  =  1000 
k  -  10-2 
=  2.34 

O  Exact  Numerical  Solution  (Ref.  5) 


c 

Fig.  2  Comparison  of  Eq.  (24)  and  Exact  Numerical  Solution 
Using  Ji,  A,  and  t  from  Ref.  5 
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Fig-  3  Degree  of  Electron  Thermo!  Nonequilibrium  in  the  Boundary  Layer  Assuming  Equilibrium  in  Free  Stream 
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A  =  1000 
k  =  10-2 

=  2.344  from  Analytical  Solution 

=  2.340  from  Numerical  Solution 

{Symbols  are  Numerical  Results 
Ref.  5) 
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C 


Fig.  4  Equilibrium  Electron  Temperature  Sheath  Profiles 
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A  =  3540 
k  =  7. 1  x  10 

D  =3.85 


=  8.74  from  Analytical  Solution 

(Symbols  are  Numerical  Results 
Ref.  8) 
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Fig.  8  Voltage  Drop  for  on  Equilibrium  and  Nonequilibrium  Sheath  and  Boundary  Layer 
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Stream  Electron  Number  Density  with  Current  Saturated 


Probe  Current 


0  2  4  6  8 

Probe  Voltage,  (Vwl  -  Vw2>  =  <*> 

TO 


Fig.  IT  Current-Voltage  Characteristics  Toro  Double  Parallel-Plate  Probe 
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Fig.  13  Variation  of  ion  Saturation  Current  with  Free-Stream  Electron  Density 


O  \ 


Ion  Schmidt  Number 
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Fig,  14  Yoriotion  of  Ion  Schmidt  Number  with  Temperature 
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APPENDIX  II 

SIMPLIFIED  SHEATH  EQUATIONS 


EQUILIBRIUM  ELECTRON  TEMPERATURE 

In  Section  III  it  was  pointed  out  that  the  integrals  in  Eqs.  (27)  and 
{28)  could  be  evaluated  in  terms  of  the  generalized  Kummer  functions. 
However,  in  the  cases  of  interest  in  this  paper  these  integrals  may  be 
simplified  and  evaluated  in  terms  of  much  simpler  functions.  Equations 
{27)  and  (28)  are: 


r0A 


a'i(t)  = 


-  t 


,t  (  *3 


3A' 


1/3 


/[/•*] 


w 


(A?'  +  1) 


e  *  dt 


(II- 1) 


oteCt> 


Yn  A 


-  E  t 


-t 


2^ 


1/3 


V  3A^ 


/  [£•■]• 


w 


Ay9  -  1 


3  e*  dt 


(II-2) 


where 


1  /  1  -  k 


i~-) 

x  1  +  k  ' 


^3 


2jx  (i  +  it) 
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Since  k  is  lCT^  or  less  for  most  conditions, 


ky0  +  1 


3(1  +  k) 


2 

ft*  -  — 


roA 


2  =  i  f  „  I 

3  '  1  +  k  '  3 


-  (4  +  Ay  )  = 
3  Z 


5  +  3k 
3(1  +  k) 


5 

3 


Using  the  above  approximations  in  Eq.  (II- 1)  gives 

J,  .  ,  ir  „  2/3 


t) 


i  (.1 

2A  '  R 


<r> 

W 


t  -  tw 


-  1 


t  ,2/3 
+  e  t 


In  Eq.  (II- 2)  we  let 

Ar2  -  4 

3 

Ayr,  -  1 


3  +  5E 
3(1  +  k) 


ftt  -  1 


which  gives 


-  2k 
3(1  +  k) 


J,  k 


w  0 


-t 


“eCt) 


r  (i)[f  (Ei(t)  -E^M 


( 1  - et*)] 


(II- 3) 


(II-  4 ) 


(II- 5 ) 


(II- 6) 

(II- 7) 

(II- 8) 

(II- 9 ) 


55 


A  E  DC-T  R-67-7  1 


where  E^(t)  is  the  exponential  integral.  Making  Eqs.  (II- 6)  and  (II-  9) 
satisfy  the  boundary  conditions 

a  -  1.0;  a  -  1.0;  C  -  1.0 

and  using  the  approximation  in  Eq.  (24)  evaluated  at  ?  =  0  gives  three 
equations  for  the  unknowns  Ji,  Rw,  and  RQ. 


If  the  asymptotic  expression  for  large  A  given  in  Eq.  (23)  is  used 
for  R0,  Eqs.  (II-6),  (II-9),  and  (24),  respectively,  reduce  to 


Jx  (1  -  E) 

2A 


(II-  10) 


and 


2E 

(3J,)273 

2A1/3  (1  +  k) l/3 


2  1  -  k  1_1 

3  1  +  E  twJ 


(11-11) 

(11-12) 


The  equilibrium  sheath  solution  for  Jj,  Rw,  and  RQ  comes  from  solving 
Eqs.  (II- 10),  (II- 11)  and  (11-12).  In  order  to  obtain  the  profiles  in  Fig.  4, 
Eqs.  ( II- 6 )  and  {II- 9 )  must  be  used. 


It  should  be  pointed  out  that  the  Rw  found  is  an  approximation,  and 
the  exact  value  of  R  at  £  =0  can  be  found  from  Eq.  (20), 


R 


w 


J^l  +  k) 
2 


+  R 


(H-13) 


CONSTANT  ELECTRON  TEMPERATURE 


The  constant  electron  temperature  case  is  treated  similar  to  the 
equilibrium  case.  Equations  (33)  and  (34)  are: 


a.(t) 


( 


IP 


) 


(11-14) 
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ox* 


(t)  =  -J1  k  t  ^3 
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Using  the  approximation 
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and  integrating  gives 
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where 


The  third  equation  is  obtained  from  Eq.  (36)  evaluated  at  ?  =  0: 

J1  (1  +  &  =■  (1  +  05o)  +  -  (  Rw2  -  Rq2  )  (11-18) 

2 

It  should  be  pointed  out  that  72  an<3  T3  are  not  the  same  as  those  in  the 
equilibrium  solution  but  are  given  in  Eq.  (30).  In  general,  approxima¬ 
tions  of  the  type  made  in  the  equilibrium  electron  temperature  solution 
cannot  be  made  in  Eqs.  (11-16)  and  (11-17).  To  obtain  the  profiles  shown 
in  Fig.  5,  Eqs.  (11-16)  and  (11-17)  have  to  be  used  in  the  form  shown. 
Solving  Eqs.  (II- 16)  and  (11-17),  evaluated  at  ?  =  1.0,  along  with  Eq.  (11-18) 
completes  the  constant  electron  temperature  sheath  solution. 
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APPENDIX  III 

BOUNDARY-LAYER  APPROXIMATIONS 


A  useful  method  of  solving  boundary-layer  equations  is  demon¬ 
strated  in  this  section.  The  Blasius  equation  is 


f*  +  fi"  =  0 


(III-l) 


where  the  superscript  prime  denotes  differentiation  with  the  variable  p. 
Equation  (III- 1 )  can  be  integrated  formally  as 


f ' 


(III-2) 


Applying  the  boundary  conditions  of  f'(o)  =  0  and  f'  =  1. 0  as  p-*®  gives 


(III- 3) 


Here  f  can  be  evaluated  for  any  assumed  f.  Since  the  exponential 
rapidly  decays  as  p  becomes  larger,  the  integration  weights  f  heavily 
near  the  wall  (p  =  0).  For  cases  of  no  slip  and  no  suction  and/or  blow¬ 
ing  at  the  wall,  the  velocity  profile  is  almost  linear  near  the  wall.  This 
implies 


f 


which  in  turn  gives 


f  = 

Substituting  this  relation 


f'  = 


c 

T1 

w 

f" 

w 

2 

2 

into  Eq.  (Ill 

1  f" 

1  w  „ 

ri 

—  *n 

3  6 

integrating  gives 


(III- 4) 
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where 


is  the  incomplete  gamma  function. 


In  order  to  check  the  accuracy  of  the  approximation,  differentiate 
Eq.  (Ill- 4 )  and  evaluate  the  derivative  at  rj  =  0 


t"  (o) 


This  gives  f^  =  0.  480  which  compares  to  f^.  =  0.  470  for  numerical 
solutions.  Hence,  Eq.  (Ill- 4)  is  seen  to  be  a  good  approximation  for 
the  velocity  profile. 


Since  f^J  is  known  for  numerical  solutions,  a  better  approximation 
to  f '  may  be  found  by  making  it  satisfy  the  condition  f  "(o)  =  f^J  =  0.  470. 
Using  this  condition  rather  than  f'  =  1.0  as  n— < gives 


f' 


(III-5) 


which  agrees  well  with  numerical  solutions  of  f '  for  all  n,  differing  by 
only  two  percent  as  rj— -  Equation  (III- 5 )  implies  that 


is  a  good  approximation  to  f"for  all  n  and  can  be  used  in  obtaining  solu¬ 
tions  for  other  quantities  in  boundary-layer  convection- diffusion  equa¬ 
tions. 
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APPENDIX  IV 

RANGE  OF  VALIDITY  OF  THE  THEORY 


An  examination  will  be  made  here  of  the  assumptions  given  in 
Section  II  to  show  the  approximate  range  of  validity  for  the  theory  pre¬ 
sented  in  this  report.  The  three  most  critical  assumptions  are  exam¬ 
ined  in  more  detail  below. 


FROZEN  FLOW 

The  flow  is  assumed  to  be  frozen  in  ionization  and  recombination 
both  while  the  particles  diffuse  through  the  boundary  layer  and  are 
convected  along  the  plate,  The  characteristic  times  for  these  two  phe¬ 
nomena  are  approximately  the  same;  therefore  we  need  consider  only 
the  convection  along  the  plate. 

The  resident  time  for  a  particle  in  the  free  stream  is 

t  =  —  sec  (IV- 1) 

CO 

The  characteristic  reaction  for  air  is  the  dissociative- recombination  of 
NO"1".  The  rate  coefficient  for  this  is  given  in  Ref.  8  with  a  resulting 
time  for  recombination  of 


T 

rec 


10  T  3,5 
e 


sec 


(IV-2) 


For  argon,  three-body  recombination  normally  dominates;  the  rate 
coefficient  from  Ref.  19  gives. 


t  =  3.01  ■  10 
rec 


12  Te 


2.94 


sec 


n 


(IV- 3) 


For  the  assumption  of  frozen  flow  to  be  valid  we  require 
T 


<  1 


rec 


which  gives  an  electron  number  density  limit  for  air 

3.5 

u  „  T 

n  <  ~  106  -e— - 

e®  L  105 


(rv-4) 
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and  for  argon. 


n  <  1.76  x  1012  /  —  T  1,47  (IV-5) 

e»  V  L  € 

In  these  equations  T  is  in  °K,  L  is  in  m,  uro  is  in  m/sec,  and  ne  is 
in  particles /m® . 

NEGLIGIBLE  CONVECTION  IN  SHEATH 


It  is  assumed  that  the  convective  contribution  to  the  current  normal 
to  the  wall  is  negligible  in  the  sheath.  If  the  convective  terms  are  re¬ 
tained  in  the  sheath,  Eqs.  (14)  and  (15)  become  (for  Te  =  T): 
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Since  the  maximum  contribution  of  convection  occurs  at  the  edge  of  the 
sheath  (£  =  1),  the  convective  terms  are  approximated  using 
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and  (see  for  example  Fig.  4) 
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The  integration  then  gives 
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The  terms  on  the  left-hand  side  of  Eqs.  (IV-8)  and  (IV-9)  are  approxi¬ 
mately  equal  to  unity.  A  comparison  of  the  convective  term  from 
Eq.  (IV-8)  with  unity  gives. 


and  from  Eq. 


Since 


(IV- 10) 


it  is  seen  that  convection  is  never  important  in  determining  the  flow  of 
the  electrons  to  the  wall.  The  ion  current  dominates  as  saturation  is 
reached  (and  the  sheath  thickens);  hence  we  shall  use  Eq.  (IV- 10)  as 
the  criterion  to  determine  at  what  point  in  the  boundary-layer  convec¬ 
tion  may  be  neglected. 


SHEATH  THICKNESS  LARGER  THAN  AN  ELECTRON-NEUTRAL  MEAN  FREE  PATH 

The  assumption  that  the  sheath  thickness  is  larger  than  an  electron- 
neutral  mean  free  path  requires  that: 
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Taking  the  mean  free  path  to  be 
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we  see  that 
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where  M  is  the  mass  of  the  neutral  particle  and  Q£n  is  the  effective 
collision  cross  section  for  momentum  transfer. 


It  is  found  that  all  the  charge  separation  is  contained  well  within 
the  sheath  when  A  >  1000.  This  fact  is  also  pointed  out  by  Chung  in 
Ref.  5.  From  the  definition  of  A  we  have 
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to  eliminate  ne^  from  the  above  equation;  hence 


(IV- 13) 


We  may  now  put  the  criteria  developed  in  the  preceding  two  sections  into 
Eq.  (IV- 13)  to  give  the  limits  on  np  . 

c  oo 

For  convection  to  be  negligible  we  find  by  substituting  Eq.  (IV- 10) 
into  Eq.  (IV- 13) 


00 


For  the  sheath  to  be  larger  than  the  electron-neutral  mean  free 
path  we  find  by  substituting  Eq.  (IV- 12)  in  Eq.  (IV- 13) 
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where  p^  is  an  atm  and  other  quantities  are  in  MKS  units. 


The  applicability  of  the  continuum  electrostatic  flat  plate  probe  to 
a  given  flow  can  be  determined  from  Eqs.  (IV-4),  (IV-5),  (IV-14)  and 
(IV- 15)  taken  collectively. 
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